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Abstract 

A search for pair production of doubly-charged Higgs bosons in the process pp — > H ++ H — > 
is performed with the D0 Run II detector at the Fermilab Tevatron. The analysis is 
based on a sample of inclusive di-muon data collected at an energy of y/s = 1.96 TeV, corresponding 
to an integrated luminosity of 113 pb _1 . In the absence of a signal, 95% confidence level mass limits 
of M(H^) > 118.4 GeV/c 2 and M(#±±) > 98.2 GeV/c 2 are set for left-handed and right-handed 
doubly-charged Higgs bosons, respectively, assuming 100% branching into muon pairs. 

PACS numbers: 14.80. Cp FERMILAB-Pub-04/045-E 
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Doubly-charged Higgs bosons appear in theories beyond the Standard Model, in partic- 
ular, in left-right symmetric models [ij, in Higgs triplet models Q] and in Little-Higgs mod- 
els 13 . The models predict dominant decay modes to like-charge lepton pairs, H — > 11. 
Pairs of doubly-charged Higgs bosons can be produced through the Drell-Yan process 
qq —>■ 7*/Z — > H ++ H . Next-to-leading order (NLO) corrections to this cross section 
have recently been calculated J]. The pair production cross sections for left-handed states 
in the mass range studied in this Letter are about a factor two larger than for the right- 
handed states due to different coupling to the Z boson. Left-handed and right-handed states 
are distinguished through their decays into left-handed or right-handed leptons. The cross 
section also depends on the hypercharge Y of the i/ ±± boson. 

The H ±i: decay width into leptons is given by Y u = (8n)^ 1 \hu\ 2 M H ±±, where is the 

n 

Yukawa coupling to leptons |2j. A limit on ft 2 „/M^±±, where is the Yukawa coupling to 
muons, can be derived from the expected contribution to the anomalous magnetic moment 
of the muon, (g — 2) M 0. This yields upper limits on the Yukawa coupling of the order 
0.1 for Mg±± = 100 GeV/c 2 . Requiring that i7 ±± bosons decay within about 1 cm of their 
production restricts the sensitivity to to approximately greater than 10~ 7 . 

Experiments at the CERN LEP collider have searched for pair production of doubly- 
charged Higgs bosons in e + e _ interactions. Mass limits for decays into muons of M(H^ i± ) > 
100.5 GeV/c 2 and M(H^) > 100.1 GeV/c 2 were obtained by the OPAL collaboration Q, 
and a limit of M(H T ,r,C) > 99.4 GeV/c 2 by the L3 collaboration Ml. Similar limits were 

n q nnn 

set for decays into electrons jo, 13] and r- leptons ja, 13, |8J . Our measurement represents the 
first i/ ±=t search at hadron colliders, and it extends significantly the range of sensitivity for 
left-handed doubly-charged Higgs bosons decaying into muons. All limits in this Letter are 
given at 95% confidence level (CL). 

The D0 Run II detector comprises a central tracking system, a liquid- argon/uranium 
calorimeter, and an iron toroid muon spectrometer Q . The central tracking system consists 
of a silicon microstrip tracker (SMT) and a central fiber tracker (CFT), both located within 
a 2 T superconducting solenoidal magnet. The SMT strips have a typical pitch of 50-80 /xm, 
and a design optimized for tracking and vertexing capability in the pseudorapidity range 
1 77 1 < 3. The system has a six-barrel longitudinal structure, each with a set of four layers 
arranged axially around the beam pipe, and interspersed with sixteen radial disks. The 
CFT has eight thin coaxial barrels, each supporting two doublets of overlapping scintillating 
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fibers of 0.835 mm diameter, one doublet being parallel to the collision axis, and the other 
alternating by ±3° relative to the axis. The calorimeters consist of a central section (CC) 
covering \rj\ up to ~ 1, and two end calorimeters (EC) extending coverage to \r)\ < 4.2, all 
housed in separate cryostats jiol ]. Scintillators between the CC and EC cryostats provide 
sampling of showers at 1.1 < \r]\ < 1.4. A muon system resides beyond the calorimetry, and 
consists of a layer of tracking detectors and scintillation counters before 1.8 T iron toroids, 
followed by two more similar layers after the toroids. Tracking at < 1 relies on 10 cm 
wide drift tubes llOl. while 1 cm mini-drift tubes are used at 1 < \r]\ < 2. 

This analysis Hj is based on inclusive di-muon data recorded between August 2002 
and June 2003. The events are triggered by requiring two muon candidates in the muon 
scintillation counters and at least one reconstructed muon using the muon wire chambers. 
The integrated luminosity is measured using two scintillator hodoscopes located on either 
side of the interaction region. For the di-muon triggers the integrated luminosity is 113 ± 
Tpb- 1 . 

Event selection proceeds in four steps. The first step (selection SI) requires at least two 
muons. Each muon used in the analysis must have a transverse momentum p? > 15 GeV/c, 
where px is measured with respect to the beam axis. The muon tracks are reconstructed 
using wire and scintillator hits in the different layers of the muon system. These must be 
combined successfully with a central track reconstructed in the SMT and CFT detectors to 
measure the muon momentum. A requirement on the timing of hits in different scintillator 
layers is used to minimize background from cosmic rays. 

The second set of selections (S2) is based on isolation criteria based on calorimeter and 
tracking information, and is designed primarily to reject background from muons originat- 
ing from semi-leptonic B hadron decays. The direction of each muon track is projected 
through the calorimeter. For at least two muons, the sum of the transverse energies of the 
calorimeter cells in an annular ring 0.1 < R < 0.4 around each muon direction is required 
to be ^ cells j E l T < 2.5 GeV, where R = \J (A0) 2 + (Ar/) 2 and <\> is the azimuthal angle. In 
addition, the sum of the transverse momenta of all tracks other than that of the muon in a 
cone of R = 0.5 around the muon track is required to satisfy ^tracks iP % T < ^-5 GeV/c. 

Selection (S3) applies to events with just two muons and requires a difference in azimuthal 
angle A0 < 0.87T. It is applied to reject Z — > fi + fi~ events and to reduce background from 
semi-leptonic B hadron decays. This selection also removes the remaining background from 
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cosmic muons. 

The final selection (S4) requires at least one pair of muons in the event to be of like-sign 
charge. These pairs are considered candidates for i? ±± — > fj, jj, decays. 

The geometric and kinematic acceptance is taken from a GEANT-based simulation 
of the detector. All other efficiencies are determined from Z — > fi + fi~ data. The single 
muon detection and reconstruction efficiencies, and the efficiency of the isolation requirement 
are measured by using one muon to tag the event and the second muon to measure the 
efficiencies. The trigger efficiency is measured by analyzing events with calorimeter-based 
triggers which are independent of the muon system. To obtain the signal and background 
rates, corrections are applied to the simulation so that it reproduces the measured efficiencies. 
The total signal efficiency for our event selection is (47.5 ± 2.5)% and does not depend on 
the mass of the doubly-charged Higgs in the mass range studied. 

TABLE I: Expected number of events for a sig nal with M(H^) = 100 GeV/c 2 , background 
events from Monte Carlo for the available integrated data luminosity, and number of observed 
events remaining after each selection cut. The simulation of Z decays includes the Drell-Yan 
contribution. 



Selection 


2 muons 


Isolation 


A4> < 0.8tt 


Like-sign 




p T > 15 GeV/c 










SI 


S2 


S3 


S4 


Signal 


9.4 


8.5 


7.5 


6.5 


Z — > /i + n~ 


4816 


4055 


359 


0.3 ±0.1 


bb 


391 


18 


3.0 


0.8 ±0.4 


Z -> T+T- 


10 


34 


2.4 


< 0.1 


il. 


3.0 


2.1 


1.5 


0.11 ±0.01 


zz 


0.1 


0.1 


0.1 


0.05 ±0.01 


wz 


0.6 


0.5 


0.4 


0.23 ±0.01 


WW 


3.5 


3.1 


1.9 


< 0.01 


Total 










background 


5254 ± 47 


4113 ±43 


368 ± 14 


1.5 ±0.4 


Data 


5168 


4133 


378 


3 
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TABLE II: Expected number of background events from Monte Carlo for the available data 
luminosity, and number of observed events remaining after each selection cut, with selection S4, 
requiring at least one like-charge muon pair, applied together with SI. Only statistical uncertainties 
are given. The contribution from WW events is negligible. 



Selection 


2 muons 


Isolation 


Acj) < 0.8vr 


(Like-sign) 


p T > 15 GeV/c 








S4 & SI 


S2 


S3 


Z — > 


0.9 ±0.3 


0.6 ±0.2 


0.3 ±0.1 


bb 


95.1 ±3.3 


4.4 ± 1.9 


0.8 ±0.4 


Z -» t+t- 
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5 
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Distributions of the di-muon mass and of A0 after selection SI are shown in Fig. ^ The 
data are compared to the sum of Monte Carlo (MC) contributions from different background 
processes. In events with more than two muons, the di-muon mass and A0 are calculated 
only for the two muons of highest px, independent of their charge. Both signal and back- 
ground events are generated with pythia 6.2 Il3f. The NNLO cross-section is used to 

n 

normalize the Z — > sample [14j . pythia does not provide a good description of the 

jet multiplicity in Z±jets events. Since the Acf> distributions are expected to be sensitive 
to the number of jets, the simulated Z — > events are re-weighted to reproduce the 

distribution of jet multiplicities observed in data. There is agreement between data and MC 
simulation, for both the normalization and shapes of the di-muon mass and A0 distributions 
(Fig.Q). The number of events remaining after each selection and the efficiencies for a signal 
of a mass of M(H^) = 100 GeV/c 2 are given in Table UJ 

The background contribution from tt and di-boson ( WZ, ZZ and WW) production is also 
estimated by MC simulation. The NLO cross section is used for ti events Higher-order 
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FIG. 1: Measured distribution of (a) di-muon mass and (b) A<p compared to the sum of MC 
background processes after selection SI. Di-muon mass and Acp are calculated for the two highest 
Pt muons, independent of their charge. 



QCD corrections to di-boson production are approximated by multiplying the LO cross 
section given in pythia by a K-f actor of 1.34. 

When the requirement of having at least one pair of like-charge muons is applied simul- 
taneously with SI, most of the background from Z boson decays is removed, and only 101 
like-sign events remain fTableITT]). Since no isolation requirement is imposed at this stage, the 
largest remaining background is from bb production, pythia is used to estimate this back- 
ground by generating inclusive jet events with a minimum transverse momentum of 30 GeV/ c 
for the hard interaction. The inclusive b quark production cross section a b (p b T > 30 GeV/ c) 
was measured by D0 to be 54 ±20 nb in the rapidity interval \y b \ < 1 at y/s = 1.8 TeV [16(. 
This cross section is extrapolated via pythia to the full y b range and to \fs = 1.96 TeV, 
and is then used to normalize the bb MC sample. 

Distributions in di-muon mass and A0 for the like-sign events are compared to the pythia 
bb simulation in Fig. El Since data and Monte Carlo are in good agreement, the efficiency 
of the A0 requirement is taken from the simulation. The data are used to determine the 
isolation efficiency Out of 101 like-sign events, five remain after applying the isolation re- 
quirement (S2). Assuming that all like-sign events originate from bb processes, the isolation 
efficiency for bb events is found to be (5 ± 2)%, and the background from bb production 



12 



in the final sample is expected to be 0.8 ± 0.3 events. Adding a systematic uncertainty of 
; measured bb cross section [16j yields a total uncertainty on the bb background 



37% 


on 


of 50%. 
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FIG. 2: (a) Di-muon mass and (b) A(f> for the two like-charge muons with the highest pt- The 101 
events remaining in the data after the selections SI and S4 (points with error bars) are compared 
to the pythia bb simulation (open histogram). The five data events remaining after the isolation 
selection are shown separately (full histogram). 

Another potential background is from Z — > decays which are not rejected by the 

A(j) < 0.87T requirement, and in which one of the muon charges is misidentified. For very 
high pt tracks, the uncertainty on the measured curvature can cause such a flip of the track 
curvature. Most of these tracks have \r)\ > 1.62, because there are fewer CFT layers in this 
region. The Z — ► simulation predicts 0.3 ±0.1 events after the final selection. We have 

also estimated the probability for charge misidentification using data. The upper limit is 
given by the ratio of like-sign (5) to opposite sign (4133) events after the selection S2 (TableslU 
and|nj), and equals 0.12%. Since 378 events remain before the like-sign requirement, then 
assuming that the charge-misidentification probability is independent of the A0 requirement, 
less than 0.5 ± 0.2 background events are expected due to charge misidentification. This is 
in good agreement with the simulation. 

The production of W bosons decaying into fiu, in association with jets, is another source 
of background, but mainly at low di-muon mass. By extrapolating to pt > 15 GeV/c the 
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FIG. 3: (a) Di-muon mass and (b) A(f> for the two like-charge muons of highest px- The data are 
compared to the sum of MC background processes after all selections. The open histogram shows 
the distributions for a left-handed, doubly-charged Higgs boson with a mass of 120 GeV/c 2 . Since 
the A(f> requirement is only applied to events with two muons, events with more than two muons 
can contribute for A(ft > 0.87T. 

steeply falling px spectrum of muons that fail the isolation requirements in di-muon events 
from a sample of W — > fiv+jets data, we estimate this contribution to be less than 0.1 
events. The expected background rate, as determined from the data, is in agreement with 
the MC simulation. 

Three candidates remain in the data after the final selection. The di-muon mass and A0 
distributions for these events are compared to the sum of MC backgrounds in Fig. EJ Two 
events have two negatively- charged muons and one positively-charged muon. Of the two, 
one has A0 = 2.48, and it has the highest like-sign di-muon mass (183 GeV/c 2 ) of the three 
candidates. The like-sign di-muon mass of the second event is 63 GeV/c 2 and the invariant 
mass of the two highest pt muons of opposite charge in this event is 91 GeV/c 2 . The third 
event has two positively-charged muons with a mass of 62 GeV/c 2 . The higher pt track may 
be interpreted of charge misidentification since it traverses the CFT in a region 

with less than 16 layers. 

Since the remaining candidate events are consistent with a background observation, i/ ±± 
mass limits are derived using the program mclimit ^1 . It provides the confidence level for 
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FIG. 4: Confidence level CLs as a function of the mass M H ±± for (a) left-handed and (b) right- 
handed doubly-charged Higgs bosons. Masses with CLs < 0.05 are excluded by this analysis. The 
mass regions excluded by LEP are also shown. The impact of systematic uncertainties is included 
in the limits. The dashed curve shows median expected CLs for no signal. 

the background to represent the data, CLr, and the confidence level for the sum of signal 
and background hypothesis, CLs+b, |l8( taking into account the expected mass distribution 
for signal and background, and the mass resolution. The mass resolution varies between 
« 7.6 GeV/c 2 for M H ±± = 80 GeV/c 2 and « 30 GeV/c 2 for M H ±± = 200 GeV/c 2 . The 
expected rate for signal as a function of the Higgs mass is determined by the NLO cross 
section J], the signal efficiencies, and the measured luminosity. The 95% CL limit for signal 
is defined as CLs = CLs+b/CLb, requiring CLs = 0.05. 

The following sources of systematic uncertainty affecting the normalization of the signal 
are taken into account: The systematic uncertainty on the luminosity is estimated to be 
6.5%. The total uncertainty on the efficiency amounts to 5%, and is dominated by the 
uncertainties on the efficiency to reconstruct an isolated muon and on the trigger efficiency. 
The uncertainty on the NLO i/ ±=t production cross section from choice of parton distribution 
functions and renormalization and factorization scales is about 10% 4]. The uncertainty on 
the background from MC is 27% (Table [TTj) . This includes the statistical uncertainty and 



the systematic uncertainty on the measured bb cross section 



16] 
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The systematic uncertainties on signal and background are taken into account in the 
limit calculation through averaging over possible values of signal and background, as given 
by their probability distributions, which are assumed to be Gaussian This procedure 
weakens the limit on the mass by about 1 GeV/c 2 . Other sources of systematic uncertainties, 
such as interpolation procedure used for the cross sections and the description of the mass 
resolution, were examined and found to be negligible. 

Figure 0] shows CLs as a function of the mass of a doubly-charged Higgs boson. The 
median expected CL$ indicates the sensitivity of the experiment for our luminosity, assuming 
that there is no signal. Taking into account systematic uncertainties, a lower mass limit of 
118.4 GeV/c 2 is obtained for a left-handed and 98.2 GeV/c 2 for a right-handed doubly- 
charged Higgs boson, assuming 100% branching into muon pairs, hypercharge Y = |2|, and 
Yukawa couplings /i w > 10~ 7 . This is the first search for doubly-charged Higgs bosons 
at hadron colliders. It significantly extends the previous mass limit [f| for a left-handed 
doubly-charged Higgs boson. 
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